Massive 2-form field and holographic ferromagnetic phase transition 
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In this paper we investigate in some detail the holographic ferromagnetic phase transition in an 
AdS, black brane background by introducing a massive 2-form field coupled to the Maxwell field 
strength in the bulk. In the two probe limits, one is to neglect the back reaction of the 2-form field 
to the background geometry and to the Maxwell field, and the other to neglect the back reaction 
of both the Maxwell field and the 2-form field, we find that the spontaneous magnetization and the 
ferromagnetic phase transition always happen when the temperature gets low enough with similar 
critical behavior. We calculate the DC resistivity in a semi-analytical method in the second probe 
limit and find it behaves as the colossal magnetic resistance effect in some materials. In the case with 
the first probe limit, we obtain the off-shell free energy of the holographic model near the critical 
temperature and compare with the Ising-like model. We also study the back reaction effect and 
find that the phase transition is always second order. In addition, we find an analytical Reissner- 
Norstr6m-like black brane solution in the Einstein-Maxwell-2-form field theory with a negative 
cosmological constant. 
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els, magnetism does not pay the central role. Yet the 
scarcity of models on magnetism is due to various tech- 
nical challenges in holographic context. In a previ- 
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izing the paramagnetism/ferromagnetism phase transi- 
tion in a dyonic Reissner-Nordatrém-AdS black brane. 
This model also was extended to realize the paramag- 
netism/antiferromangnetism phase transition by intro- 
ducing two real antisymmetric tensor fields with inter- 
action between them and they are coupled to the back- 
ground gauge field strength [24]. In Ref. [25], by combin- 
ing the holographic p-wave superconductor model [26] 
and the holographic ferromagnetism model, we studied 
the coexistence and competition of ferromagnetism and 
p-wave superconductivity. It is found that the results 
depend on the self-interaction of magnetic moment of 
the complex vector field and which phase (superconduc- 
tivity or ferromagnetism) appears first. We noted that 
the ferromagnetic superconductivity has been also dis- 
cussed in Ref. [27] by introducing two SU(2) gauge fields 
in the bulk. In that model superconductivity and fer- 
romagnetism happens simultaneously and the magnetic 
susceptibility is finite at the ferromagnetic critical tem- 
perature. However, the ferromagnetic p-wave supercon- 
ductor in heavy fermion systems such as UCoGe, URhGe 
and UGeg, shows the two critical temperatures are dif- 
ferent in general [28] and magnetic susceptibility diverges 
at ferromagnetic critical temperature [29]. 

The reasons of using an antisymmetric tensor field 
to model ferromagnetic phase transition in holographic 
setup are as follows. On the viewpoint of symmetry 
breaking, ferromagnetic phase transition breaks the time 
reversal symmetry spontaneously in low temperature (if 
spatial dimension is more than 2, it also breaks spatial 
rotation symmetry), which in general is not associated 
with other symmetry breakings such as U(1), SU(2) and 
so on. Thus the dual operator does not carry U(1) or 
SU(2) charge, which implies that we need to use a real 
field corresponding to such an operator. From the point 
of view of covariance, magnetic field is the spatial compo- 
nent of a SO(1,3) tensor Fu, so magnetic moment should 
be also a component of a tensor. Furthermore, according 
to the origin of magnetic moment, magnetism of mate- 
rial is controlled by the intrinsic magnetism of electrons 
which is obtained from the Lagrangian for electrons, 


Le = (bd — mb) — ej" A, — 5 bot bE wy: (1) 


where A, is the gauge potential, F,, is the gauge 
strength field and j” is current density which is defined 


> £ 
-zz Y( ð” — 0"). Then the magnetic moment 


density N is the spatial component of an antisymmetric 
tensor field, 


as j” = 


N? = êj MIE, MH” = — (oy). (2) 
4m 
Here i,j,k = 1,2,3. Thus we see that an effective field 
to describe magnetic moment in the boundary field in a 
covariant manner needs an antisymmetric tensor and its 
spatial component corresponds to the magnetic moment. 
These considerations can be generalized into 2+1 dimen- 
sions, where magnetic moment should still be regarded as 


the spatial components of an antisymmetric tensor field. 
Then according to the spirit of AdS/CFT, we can use 
an antisymmetric tensor field in the bulk to dual such a 
tensor operator. 

Recently, it has been found in Ref. [30] that the origi- 
nal model in [23] contains a ghost and a new model has 
been proposed, there the exterior differential is used to 
replace the covariant derivative in the kinetic term of the 
antisymmetric tensor field. With such a simple modifica- 
tion it has been shown that there does not exist any ghost 
and causality violation does not appear in the new model, 
while the key results in the original model are kept qual- 
itatively. While in Ref. [30] we paid our main attention 
to the health of the model and showed that the sponta- 
neous magnetization can happen when the temperature 
gets low enough, in this paper, we are going to study 
this model in some detail with some directions. Before 
going on, we mention here that in this paper we assume 
the physical phase is homogeneous and will not consider 
inhomogeneous phases. This assumption is just for sim- 
plifying our discussions in technology. The possibility 
whether the inhomogeneous phase could appear sponta- 
neously in this model under the homogeneous boundary 
conditions is left to study in the future. 

This paper is organized as follows. In section II, 
we will describe the holographic model. We will give 
our ansatz for matter field and derive equations of mo- 
tion in section III. In this section, we will show that 
there exists an analytical black hole solution different 
from the AdS Reissner-Nordtrs6m solution in this model 
and discuss associated properties of the solution. In 
sections IV and V, we will investigate the paramag- 
netism/ferromagnetism phase transition in two different 
probe limits. One is to neglect the back reaction of 
the 2-form field to the black brane geometry and to the 
Maxwell field, and the other to neglect the back reaction 
of both the Maxwell field and the form field. In the for- 
mer case, we will also calculate the off-shell free energy of 
the holographic model near the phase transition point by 
using the Sturm-Liouville eigenvalue method and obtain 
a free energy form like the Ginzburg-Landau one, and 
make a comparison with the Ising like university class 
model. In the latter case, we calculate the DC resistiv- 
ity in the ferromagnetic phase, which shows the behavior 
of the colossal magnetic resistance effect in condensed 
matter physics. We study the full back reaction effect in 
section VI by solving the full equations of motion of the 
model and find that the phase transition is aways second 
order. The summary and some discussions are included 
in section VII. 


Il. THE MODEL 


In this paper, the model we are considering is just 
general relativity with a negative cosmological constant 
A = —3/L?, a U(1) field A, and a massive 2-form field 
M,, in 4-dimension space-time. The ghost free action 


reads [30] 
1 4 2 
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where L is the AdS radius which is set to be unity and 
2k? = 167G is associated to the gravitational constant 
in the bulk, which will be set to be unity in the follow- 
ing. g is the determinant of the bulk metric guv. dM 
is the exterior differential of 2-form field M,,. m? 4 0 
is the squared mass of 2-form field M,,, and must be 
greater than zero, which will be explained shortly. À and 
J are two real model parameters with J < 0 in order the 
magnetization to happen spontaneously. \? characterizes 
the back reaction of the 2-form field M,, to the back- 
ground geometry and to the Maxwell field strenghth 1. 
V(M,„v) is a nonlinear potential of the 2-form field. It 
describes the self-interaction of the polarization tensor, 
which should be expanded as the even power of M,,,. In 
this model, we take the following form, 


V(Myv) = “M, Me’)? = [*(M AM). (4) 
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Here * is the Hodge-star operator. The choice of nonlin- 
ear potential is not unique. We choose this form just for 
simplicity. As shown in Ref. [30], this potential shows a 
global minimum at some nonzero value of p. 

By varying action (3), we can get the equations of mo- 
tion for the matter fields and gravitational field as 


V7 (dM)ruv — M My — J Mrs M7?) (* My) = Fu, 


A2 
V” (Fav + go) =0, (5) 


1 3 
Raw = gPIw = T29. = Tyv- 
The energy-momentum tensor T),,, reads 
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In the AdS/CFT correspondence, a hairy black hole 
with appropriate boundary conditions can be explained 
as a condensed phase of the dual field theory, while a 
black hole without hair is dual to a normal phase. Clearly 
when M,,, vanishes, the model admits the AdS Reissner- 
Nordström (RN) black brane solution, which corresponds 


1 Note that here \ can be also understood as the measurement of 
the coupling between the tensor field My, and the Maxwell field 
strength Fuv by rescaling the tensor field and the parameter J. 


to the normal phase in the dual field theory. When M,,, 
appears, we will see that there exists an analytical black 
brane solution differing from the AdS-RN solution. The 
new solution corresponds to the normal phase rather than 
the usual AdS RN solution. When we lower the tempera- 
ture, the system exhibits an instability which triggers to 
break time reversal symmetry spontaneously as well as 
spatial rotation symmetry since the condensate will pick 
out one direction as special (if spatial dimension is more 
than 2) and the paramagnetism/ferromagnetism phase 
transition happens. 


II. ANSATZ AND TRIVIAL SOLUTION 


As we will discuss the full solution to the action (3) 
including the back reaction to the spacetime geometry, 
we start with the following ansatz for the metric 


dr? 
r? f(r) 


ds? = —r? f (rje dt? + + r°(dz? + dy’), (7) 


and take the self-consistent ansatz for polarization field 
and U(1) field as 


Mav = —p(r)dt A dr + p(r)dz ^ dy, 
A, = o(r)dt + Bady, (8) 


with some real functions f(r), a(r), (r), p(r),and p(r). 
The bulk field B is a constant magnetic field, which can 
be regarded as external magnetic field in the dual bound- 
ary theory. We will denote the position of the horizon as 
rp, and the conformal boundary will be at r > oo. Since 
we would like to study a dual theory with finite chem- 
ical potential or charge density accompanied by a U(1) 
symmetry, we turn on A; in the bulk. 

Suppose the black brane horizon is at r, with f(ra) = 
0. According to gauge/gravity duality, the Hawking tem- 
perature of the black brane is identified with the temper- 
ature of boundary thermal state, which is given by 


(F(Y e 
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T= 


and the thermal entropy S is given by the Bekenstein- 
Hawking entropy of the black brane 


S = 4TA = 4rr? Va, (10) 


where A denotes the area of the horizon and V2 = f dzdy. 
Put the ansatz (7) and (8) into equations (5), the in- 


dependent equations of motion read 
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where a prime denotes the derivative with respect to r. 
The first three equations are for the polarization field 
and Maxwell field, and the last two are the two indepen- 
dent components of gravitational field equations. In fact, 
there are three nonzero components of gravitational field 
equations, but only two of them are independent due to 
the Bianchi identity. 

We are interested in the black brane configurations 
which have a regular event horizon located at ra. There- 
fore, in addition to f (ra) = 0, one must require (r) = 0 
in order for g#” A,A, being finite at the horizon. We re- 
quire the regularity conditions at the horizon r = rp, 
which means that all functions will have finite values at 
rp, and admit a series expansion in terms of (r — rp). 
Then, at the horizon, we have the following relations, 


AJ p? 
o=0, p = (m° — ri )p, 
p(m? + 4Jp?e~*) + B 
= F i (12) 
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The black brane solution should be asymptotically AdS. 
Thus one has the following asymptotic solutions near the 
AdS boundary, 


H3 + B?. 
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where 6 = /1+4m? and the dots stand for the higher 
order terms of 1/r, p+ and fo are all constants. The 
Breitenlohner-Freedman (BF) bound requires m? > —4 
according to the asymptotic solution of p. According to 
the AdS/CFT dictionary, up to a normalization, the co- 
efficients u and o are directly related to the chemical 
potential and charge density in the dual system, respec- 
tively. The magnetic moment density in the dual theory 


is defined by the integration [23, 30] 


X2 oo pel? 
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N= dr. (14) 


Note that there is an additional restriction on m? when 
we treat the magnetic field B as the source of p. In order 
to keep p damping, we have to impose boundary con- 
dition for p such that p} = 0 and a restriction on the 
parameter such that m? > 0 [30], otherwise the integra- 
tion (14) will diverge or the equation for p can not be 
linearized near the AdS boundary. These two restric- 
tions can be understood in another way: when B = 0, 
the constant p} should be viewed as the source of the 
corresponding operator in the boundary field theory, ac- 
cording to AdS/CFT duality. In order the symmetry to 
be broken spontaneously, one has to turn off the source 
term. 


In addition, in order to make p condense in the case 
without external magnetic field when the temperature is 
low enough, the model parameters should violate the BF 
bound for p near the horizon, whose geometry is an AdS2 
for an extremal black brane. 

To see how this requirement restricts the parameters, 


we first consider the solution for (11) in the case of p(r) = 
B = 0 and a(r) = V(M) = 0. The equations (9) become, 


mp ~ o = 0, 
2 , 
o" + “9 xE + =) = 0, (15) 
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The equations admit an analytical solution when m? Æ 0, 
i.e., 


o(r)=H(L—ra/r), p) =F, (16) 
f(r) =1 cee + Eri, p= ar = a 
Then the temperature of the black brane is 
= =(3 - ñ). (17) 


Compared with the usual planar AdS RN black brane 
solution [31], we see that u? in the metric of usual AdS 
RN solution is just replaced by ñ in the new solution. 
Namely, the new solution has the same form as the planar 
AdS RN solution with the same Maxwell field ¢(r), but 
also with a nontrivial profile of p(r) for the massive 2- 
form field. 


Then we can easily calculate the free energy Q and the 


charge density o of the system as ?, 


Q = Q(T, p) = —2Vər}(1 + F) 


2 (09 A2 (18) 
o =o(T, u) = v; (5). = 2(1 Td hte 

We see that the properties of this black hole solution de- 
pends on the value of 1 — A?/4m?. When 1 — A?/4m? = 
0, the solution is just the planar AdS Schwarzschild 
black brane with nonzero U(1) gauge field and polar- 
ization field. To investigate the physical properties when 
1— à? /4m? Æ 0, it is useful to compute the partial deriva- 
tive of charge density with respect to chemical potential, 


Oo 1 (A?Q X 1+5 
Ge V2 (Sr), 780 pads 09) 


It is easy to see that Eq. (19) is positive when 1 — 
à? /4m? > 0. However, when 1—A?/4m? < 0, the value of 
(19) can become negative, which leads the dual boundary 
system to be chemical instability. To understand it, one 
can image that a box is submerged into the environment 
with fixed temperature T and chemical potential jw in 
equilibrium state. The energy and charge can exchange 
between the interior of the box and environment through 
the wall (see Fig. 1). Now because of thermal fluctuation, 
the chemical potential in the box then is u — du < p, 
which leads to some charge coming into the box. How- 
ever, because Eqs. (19) is negative, adding the charge 
of the system will decrease the chemical potential of the 
system. Then we see that the chemical potential will de- 
crease again, which will lead that more charge come into 
the box. So the charge inside the box will be more and 
more and the system is unstable. From Eq. (19), we can 
find the dual boundary system is chemical stable region 
only when temperature and chemical potential satisfy’, 


V27 m 
5 (2 -—1)T 


4m? 


<7. (20) 


If Eq. (19) is positive, we see that adding charge will 
increase the chemical potential, so after some time, the 
system can be in equilibrium again. In fact, in the case 
of 1 — \?/4m? > 0, the solution shares all the prop- 
erties of usual AdS RN black brane except we have to 
replace u? by ji. But in the case of 1 — A?/4m? < 0, 
the properties are very different from those of usual AdS 
RN black brane. For example, the temperature gets in- 
creased when we increase the chemical potential while 


2 The total charge density ø can be also calculated from the second 
equation in (5), which gives the same value as the one in (18). 
This also confirms that the chemical potential of the system is 
given by p. 

3 It is worth noting that this condition is only valid in grand canon- 
ical ensemble. In other ensembles, the stable condition is differ- 
ent in general. 


FIG. 1. The schematic diagram for chemical instability. 


fixing the horizon radius and there does not exist zero 
temperature entropy. 

In this paper, we only focus on the situation of 1 — 
d?/4m? > 0, which can give a chemical stable dual 
boundary system and zero temperature black hole so- 
lution whose IR geometry has the form of asymptotic 
AdS geometry. With this, we consider the possibility 
of spontaneous symmetry breaking of this system in low 
temperature. This can be analyzed by solving the follow- 
ing equation in the background (16), 


p" + A joe 4Jp? 
F CAF T maf 


The condition of instability for p in some low tempera- 
ture is that the BF bound in AdS; is retained but the 
BF bound in AdS is violated. Near the horizon for an 
extremal black brane, the geometry is asymptotic AdS2. 
In this region, the asymptotic solution for equation (15) 
is, 


jp=0. (21) 


p = pyr AV ita) /2 + pp Ovi) /2 
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So the conditions for spontaneous condensate are, 
m >0, m?(4m? -— X) > 0, (23) 
Pp m? 8J 1 
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6  m?(4m?-— X2 
One can immediately see that these inequalities have so- 
lution only when J < 0. 

Under the restriction of (23), there is a critical temper- 
ature, lower than which the nonzero p begins to appear. 
Then the full coupled equations of motion do not admit 
an analytical solution. Therefore, we have to solve them 
numerically. We will use shooting method to solve equa- 
tions (11). In order to find the solutions for all the five 
functions F = {p, ¢,p,a, f}, where we can directly solve 
(r) = ie dr(m? — 4Jp?/r*)p and then put it into the 
equation of @” and obtain an equation for p, we must 
impose suitable boundary conditions at both the horizon 
r = rp, and conformal boundary r > oo. There are two 
kinds of scaling symmetries which are useful when we 
perform numerical computations: 


ror, (t,z,y)—> (t,2,9)/a, (24) 
(p, B) > 07(0,B), $> a, T> aT. (25) 


and 
a> atà, e+ ere", +e-3t, p> e?g. (26) 


Under the above two scaling symmetries, we finally have 
four independent parameters {a(rn), p(ra), p(ra)} and 
horizon radius r, at hand. In this paper, we will fix the 
boundary chemical potential to be unitary. When these 
four parameters are given, we can integrate the equations 
out of the horizon to obtain the whole solutions. When 
we perform the numerical computations, we can first set 
{ra = 1,a(rn) = 0}. Then for a given p(r;,), we can ad- 
just the value of p(r;,) to match the boundary condition 
p+ =0 at r— oo. After solving the coupled differential 
equations, one should use the second scaling symmetry 
(26) to satisfy the asymptotic conditions a(co) = 0. Then 
we use the first transformation (24) to fix the chemical 
potential for each solution the same. By this method, we 
get one-parameter solution p(r;T) for T < Te. And then 
we can get the behavior of magnetic moment density N 
with respect to temperature T. 

When the parameters satisfy the inequalities (23), the 
solution of p Æ 0 will appear even in the case of B = 0. 
Because under the time reversal transformation, B —> 
— B, in order to make action be invariant, we have to have 
p — —p. So when B = 0 but p Æ 0 in the source free case, 
the time reversal symmetry is broken spontaneously. 

In order to see the main properties of the model, let 
us first study the model in probe limit for simplicity. 
Here we may consider two kinds of probe limit. The first 
one is to take the model parameter A — 0 as we did in 
Ref. [23]. This case is just to neglect the back reaction of 
the massive 2-form field to the black brane geometry and 
the Maxwell field, but considers the effect of the Maxwell 
field to the background geometry. This probe limit cor- 
responds to the case in such materials that their elec- 
tromagnetic response properties are very weak compared 
with the external field and have little effects on the their 
transport properties. In other words, in this probe limit, 
the influence of external field on the materials is consid- 
ered, but we neglect the back reaction of electromagnetic 
response to the external field and structures of materi- 
als such as crystal structure or energy band. The other 
is to neglect all back reaction of matter fields including 
the Maxwell field to the background geometry. In this 
probe limit, the interaction between the electromagnetic 
response and external field is taken into account so that 
we can study how spontaneous magnetization influences 
the electric transport, but they both have little influence 
on the structures of materials. We will study two cases 
in the following sections separately. 


IV. PROBE LIMIT IN THE CASE OF \—- 0 
A. Spontaneous magnetization and susceptibility 


Let us first investigate the spontaneous magnetization 
in the limit of A > 0. In this limit, we neglect the back 


reaction of polarization field to the gauge field and back- 
ground geometry. The background geometry and the 
Maxwell field can be taken as 


=al- i) fr) =1- FF +4, 
T=7(3-1), (27) 


where we have set the horizon radius ry, = 1. In this case, 
the equations of the polarization field read 


TR i m? +4Jp’? = B 
ae ae Pp 
4 Jp? lu 
(m? — a )p— a= 0. (28) 


Note that the expression of p(r) can be solved directly. 

We put it into the equation of p(r), and get 
p+ f p = ane lp 5 
f r2f (m2r2 — 4g \2p2 f r2f 


= 0. (29) 


This equation shows that p will be spontaneously con- 
densed below a critical temperature only when m? > 0. 
When m? > 0, near the critical temperature where p is 
very small and therefore the term p? can be neglected 
in (29), thus increasing the chemical potential* will de- 
crease the effective mass at the horizon (note that J < 0), 
which leads that can be condensed spontaneously below 
a critical temperature. ° Furthermore, the restrictions on 
the parameters in the case of A > 0 are, 


2 
Lo m 2J 1 2 

= —+4+— <-- > 0. 30 

k 6 7 m4 = 4’ < po 

In what follows, since different parameters will give sim- 

ilar results, we will choose m? = —J = 1/8 as a typical 


example when we perform numerical computations. 

To study the spontaneous magnetization, we need set 
B = 0. Since near the critical temperature, p is very 
small, the nonlinear term of p can be neglected in (29).To 
find the critical temperature, we can solve the linearized 
equation of p with initial condition 


,_ MË +4J(3 — 47T) 
p= 4rTm4 


at the horizon r = 1 and boundary condition p4 = 0 at 
r — oo. Without loss the generality, we take the initial 
value of p at horizon to be unity, and treat T as the 
shooting parameter to match the source free condition. 
There will be many solutions for shooting parameters T, 
we choose the highest one as the critical temperature T}. 


(31) 


4 This implies that the temperature is decreased in grand canonical 
ensemble [see (27)]. 

5 If m? <0, p will get spontaneously condensed at a temperature 
higher than a critical temperature. This solution is unstable and 
we will not discuss it any more. 
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FIG. 2. Top panel: The value of p(r;,) in different tempera- 
ture. Bottom panel: The magnetic moment N as a function 
of temperature. Here we choose the model parameters as 
m? = —J = 1/8. The critical temperature T./p ~ 1.78. 


We find the critical temperature Te /u ~ 1.78 for the case 
of given model parameters. 

When the temperature is lower than the critical tem- 
perature Te, in order to examine whether the polariza- 
tion field p can make spontaneous magnetization when 
the external magnetic field B = 0. We plot the value of 
p+ versus p(rp;,) at the horizon in the top panel of Fig. 2. 
Each curve in this plot corresponds to different tempera- 
ture. This plot shows a typical example in the high and 
low temperature cases which correspond to the red and 
the black lines, respectively, while the case with the crit- 
ical temperature is shown by the blue curve. In the case 
of high temperature, we see that the curve has no inter- 
secting point with horizontal axis except a trivial point 
at the origin, which corresponds to a trivial solution with 
p = 0 and it describes the paramagnetic phase. On the 
other hand, when the temperature is low enough, we find 
that there exists a nontrivial solution which locates at 
p(rnz) # 0. This solution breaks the time reversal sym- 
metry and the system gets into a ferromagnetic phase. 
As a result, we see that the model indeed can give rise 
to a paramagnetism/ferromagnetism phase transition in 
the case without external magnetic field. 

Furthermore, when the temperature is lower than the 
critical temperature T., we have to solve the equation 
(29) with the condition of p+ = 0 by shooting method. 
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FIG. 3. The behavior of the inverse susceptibility density in 
the paramagnetic phase near the critical temperature when 


m? = —J = 1/8. Here we set 2%? = 1 for convenience. 


The spontaneous magnetization N is defined by the inte- 
gration (14) with a = 0. The spontaneous magnetic mo- 
ment with respect to temperature is shown in the bottom 
panel of Fig. 2. In addition, since the magnetic moment 
of the polarization field obtains an expectation value, the 
time reversal symmetry is broken spontaneously. 

By fitting this curve near the critical temperature, we 
find that there is a square root behavior for the magnetic 
moment versus temperature, which is a typical behavior 
for a second order phase transition, specifically, for m? = 
—J = 1/8, we have 


N? JA u? ~ 4.910(1 — T/Te). (32) 


This gives the critical exponent 1/2, the same as the one 
from the mean field theory. 

Except for the magnetic moments which is one of the 
characteristic properties of ferromagnetic material, an- 
other remarkable one is the behavior of susceptibility 
density of the material in the external magnetic field. 
The static susceptibility density is defined by 


_ ON 

x= Jim, oe (33) 
When we turn on the external magnetic field B, the func- 
tion p is nonzero in any temperature. In order to compute 
the susceptibility density, we need to shoot for the bound- 
ary conditions with one parameter p(r;,) for equation (29) 
under the given external magnetic B and temperature T. 
From the definition of x, which involves only the behav- 
ior of B > 0, so the computation can be simplified in the 
following way. 

In the case of T > T,, considering the result p = 0 
when B = 0, we except that p is in the same order as 
B, then p° is as the same order as B?. So in the case 
with weak magnetic field, we can neglect nonlinear terms 
in equation (29). In that case, we just need to solve the 
“linearized” equation 


"y A m? 4J? B 


phe ua a a 


Because p x B, we can set B = 1, Solving the equation 
(34) with the boundary condition p+} = 0, the suscepti- 
bility density then can be computed by 


ne a x f spar. (35) 


In Fig. 3, we show the magnetic susceptibility as a func- 
tion of temperature and find it satisfies the Curie-Weiss 
law of ferromagnetism in the region of T > T+. Con- 
cretely, for the chosen model parameters, we have 


A? /ux œ 4.0499(T/T, — 1). (36) 


So we can conclude that the dual system is in a paramag- 
netic phase in high temperature and ferromagnetic phase 
in low temperature. The model can describe a param- 
agnetism/ferromagnetism phase transition. We will see 
later that the exponents in Eqs.(32) and (36) are exact 
and can be obtained by analytical methods. 


B. Holographic Ginzburg-Landau formulation 


In action (3), two parameters appear in the La- 
grangian. One is the mass of the polarization field, which 
corresponds to the conformal dimension of the dual op- 
erator by the standard dictionary. The other is the self- 
interaction coupling constant J. So a natural question is 
what the meaning of this constant is in the dual boundary 
theory? In this subsection, we will try to give a quanti- 
tative interpretation. Since we are working in the probe 
limit, the geometry and external Maxwell field are fixed, 
this leads to a simplification to compute the partition 
function of the bulk theory. This is the reason that we 
consider this kind of probe limit here. Our method is 
to compute the effective grand thermodynamic potential 
in both sides and to equate them (which are equivalent 
to compute the partition functions in both sides and to 
equate them). Then we can “read off” the meaning of J 
in the dual boundary theory. 

We assume that the the spontaneous magnetization in 
the boundary theory can be described by so-called Ising 
universality class model at least in the vicinity of critical 
temperature, whose Hamiltonian is, 


1 
=I PF) FO) + OKI IP, 
rr r 

(37) 
where r,r’ are the positions of lattices, J (r — r’) > 0 
is the exchange integration of lattices r and r’, s(r) the 
z-component of spin at lattice r, K is a positive constant 
which corresponds to the deviation from the Ising model. 
When K —> oo, Hamiltonian (37) reduces to that of the 
usual Ising model. 

In this paper, the dual boundary is a (2+1)- 
dimensional spacetime, i.e., a film system. Then the spin 
(or magnetic moment) in fact is just a pseudoscalar. So 
we have s = s. To compare with our holographic model, 


it is convenient to change Hamiltonian (37) into the one 
in continuous limit, 


H =| ea [5 FPR s(0) ¥o(0) (38) 


(J +2K)s(x)? + Ks(x)"] , 


where | is the lattice spacing and J = >>, J(r), RZ = 
$r? I(r). The summations for J and Rẹ, are only in 
one crystal lattice. We see that the Ising like Hamiltonian 
is a \d* theory if we use 3 = sy J Rg. From Appendix A, 
we can find that the grand thermodynamic potential in 
the mean field approximation in the high temperature 
limit reads, 


a 1 As YET ee) Aneel 
wi _ c 39 
tana aa TF (39) 
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Here we have used 5,; to represent the classical value 
of s. We see that there is a critical temperature T, = 
4rm?/Asyg. When T > T., the thermodynamic equilib- 
rium phase corresponds to Sea = 0, while T < Te, the 
thermodynamic equilibrium phase has Sa, Æ 0. 

Now let us consider the gravity side of the holographic 
model to compute the grand thermodynamic potential. 
It is convenient to make a coordinate transformation by 
z = rp/r. As p can be solved directly, then we put it into 
the equation of p(r), and get 


4J p22? i= B 
(m2 — 4J p2z4)2 f p= 2p 
(m? — 4. p?z*)p — uz? = 0. 
(41) 


As we will care about the behavior of T —> Te, the value 
of p will be a small quantity near the transition point. 
In this case, we can make a Taylor’s expansion on the 
nonlinear term of p in Eq. (41) as, 


AT p22? O AT p22? A 32J? u? p72 
(m2 — 4J p?z4)2 T mt ' mê 


+ O(p") (42) 


Neglecting the high order terms, Eq. (41) can be rewrit- 
ten as 


(43) 


Jy = —32J? u? /m® < 0. 


Up to the order of p*, the part of polarization field in 
action (3) can be written as, 


Zh 


S(T, B; p) 
XV 


z? 1 2 
= (z fee + zfe") 


-f dz 
0 


which is a function of T and B, but a functional of p. 
The asymptotic solution for (43) is 


° (44) 


PF Jf g a4 
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p=p-—,, with p= ps2 (1+8)/2 + g_z7078)/2, (45) 


The source free condition is p} = 0 as z + 0+. Under 
this, the grand thermodynamic potential or free energy 
in grand canonical ensemble Q is, 


Q(T, B; p) = Q(T, B; p) V2 


Zh 
= xv: | dz 
0 


According to thermodynamic relationship, 
dQ(T, B) = —SdT — V2NdB 


A oo P (47) 


(46) 
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It seems that the magnetic moment should be, 


2 à a 
PE. (mn =x f pdz. (48) 
T,p 0 


Vo OB 


However, in our previous papers, we defined the magnetic 
moment as, 


N/? = — f Pdz. (49) 
0 2 


There is a difference factor 1/2 between (49) and (48). 
Now we show that the expression (48) is not true. The 
reason is as follows. The relation (47) is on-shell, while 
the equation (46) is off-shell. In order to obtain the differ- 
ential relation of Q(T, B; p) with respect to (T, B, p), we 
need to use the Euler homogenous function theorem. We 
should first note that under the scaling transformation 
z— kz, (t,x, y) > k(t,x,y), we have Q(kT, k? B; k?p) = 
kQ(T, B; p), which gives, 


OQ OQ 
p,B pT 
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Submitting (46) into (50) and considering the on-shell 
condition 6Q/dp = 0, we find that, 


1 (dX(T a zh p 
nix = -y (A =- f =dz. (51 
/ V2 3B T,on-shell 0 2 E ( ) 


(50) 


This is just the definition (49) . 

The key step for computing the grand thermodynamic 
potential is to structure the Sturm-Liouville problem®, 
which is the following ODE: 


Lon 

G) 

= 5 {-E| POE] +a) on = anon 
(52) 


II 


Pop, 
W 


with the boundary conditions: 
(a) At z = zn, f(zn) = 0, |Pn(zn)| is required to be finite; 
(b) At z + 0*, we impose p,,(0) = 0. 


The weight function w(z) can be an arbitrary positive 
continuous function in the region of [0, zn]. From a prac- 
tical point of view, we choose weight function such that 
the values of A, will not influence the asymptotic be- 
haviors of equation (52) when z > 0*. There are many 
choices for weight function. Here we choose w(r) = 2" 
with an integer k > 2. 

Once note that the asymptotic solution when r — co 
for equation (52) is, 


Pn = pee + p-27 0-8/2 (53) 


one can find that the boundary condition (b) corresponds 
to p} = 0. Let L?([0, zn], w(2), dz) be the Hilbert space 
of square integrable functions on [0, za], i.e., 


L?([0, zn], w(z), dz) 


f w(z)|h(z)|?dz < co} (54) 
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with the inner product 
(hy, he) = f w(z)hy(z)ho(z)dz, (55) 
0 


and D be the subspace of £*((0, zn], w(z), dz) that satis- 
fies the boundary conditions of (a) and (b), i.e., 


D = {Vh € L?((0, zn], w(z), dz) |h € C?(0, zn] , (56) 
h(0) =0, |h(zn)| < oo}. 


Then we can prove that P is the self-adjoint operator on 
D, i.e., 


Yhı, ho € D, (hı, Pho) = (Phi, ho). (57) 


6 The method is similar to the one used in Ref. [32], but is 
completely different from the Sturm-Liouville (SL) eigenvalue 
method in Ref. [33], there the precision depends on the trial 
function one chooses. 
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FIG. 4. The function y(A) = arctan[10p+(A)] with re- 
spect to A at temperature of T œ 2.007. (left one) and 
T ~ 0.4917.(right one) in the case of m? = 1/8. The x- 
coordinate values of zero points give the eigenvalues An. Here 
Te ~ 1.7766p. 


According to the properties of SL problem, the solutions 
of (52) form a function basis on D with which one can 
expand any functions belonging to D, i.e., 


(Pn; Pk) = Onk; (58) 


and 


Yh € D, Hen} C R, h(2) = X` cnpn(z) (59) 


with cn = (Pn, h). 

Note that all the parameters in equation (52) depend 
on temperature, so do the eigenvalues An and eigenfunc- 
tions pn. Indeed, the minimal eigenvalue, i.e., the first 
eigenvalue À; is the function of temperature. There is 
a critical temperature Tę, at which we have A; = 0. 
We can find that when T > Te, 41 > 0. In order to 
show this, we take m? = 1/8 as an example to plot 
p+ with respect to A. The boundary condition (b) is 
equivalent to p4(An) = 0. Fig. 4 shows the value of 
y(A) = arctan[10p4(A)] with respect to A1/? (it is just 
for convenience). The x-coordinate values of zero points 
(y = 0) give the eigenvalue An. We can see that all the 
eigenvalues are positive when T > Tẹ and A, < 0 when 
T < Te. For any temperature, the equation has infi- 
nite eigenvalues and the smallest eigenvalue exists, which 
shows the system is stable. 

Let us now turn our attention to the free energy (46). 
For convenience, we will use scaling transformation to set 
Zp = 1 in the process of computation, and then transform 
into the case of fixing chemical potential in the final re- 
sults. 

Let p(r) = p(r) + B/m? be any function configuration 
belonging to D, in which p(r) dose’t need to be the so- 
lution of EoM (43). We can use the eigenfunction pn to 
expand p(r) and magnetic moment as, 


co co B 
p= 2 enPn oe p= 2 CnP = 
n= n= (60) 
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where cn and Nn are coefficients, defined as, 


1 1 
Cn =} wppndz, Ny =|} Pndz. (61) 
0 0 


Here we have assumed that {pn} is an unit base. Then 
the variational principle of Q(T, B;p) underlying the 
equations of motion, or finding a solution of EoM (43), 
is equivalent to minimizing Q(T, B,c,) with respect to 
Cn S. 

Let us consider the case of spontaneous magnetization, 
i.e., the case with B = 0. In this case, we have 


1 
Q(T, cn) = » f dz [PP Fy28 o'/4] , (62) 
0 


with p = J>] Cnn. Using the orthogonal relationship, 
we have, 


If T > Te, then A, > 0. Because of J; < 0, we can find 
that Q(T, cn) > 0. The minimization of Q(T, cn) = 0 is 
achieved only when cn = 0, i.e., p = 0. So the nonzero 
solution appears only when A; < 0, i.e., T < Te. This is 
just what we have obtained in the pervious section. 

When T > T; , we can set Ay = ao(T/Te—1) with ap > 
0 and assume that the off-shell solution is dominated by 
the first term in (60) only, i.e., [ci] > cn for n > 2 in 
(61). As a result, we have, 


~ 1 Jyh f' 
A? O(T, cn) & =Are? — a | dzp;2°, 
0 


5 i (64) 
~ zo(T/T. - 1) — Tpctay 
with a, = if p$z8dz|r=r. > 0 and, 
N ~ =X c N1 /2. (65) 


Put (65) into (64), we can obtain, 


A(T, cn) ~ AT, N) 


T 2ao 2 _ 

~ 2N? (T/T. —1)N* + N 
We can see that this is just the Ginzburg-Landau (GL) 
theory of ferromagnetic model. 

Now the grand thermodynamic potentials of the holo- 
graphic model in the gravity side and in the boundary 
theory side are in hand, we can use the AdS/CFT dual- 
ity to relate them, which gives, 


Qorr = Qgravity- (67) 


Here we should mention that the dual theory in the 
boundary being an Ising-like model is an assumption. 
With this equality, At the critical temperature, compar- 
ing (66) with (39), we can obtain, 


— N’ w ASi (68) 


Note that in materials, the spontaneous magnetization 
is proportional to the expectation value of z-component 
of spin, i.e., N x Sa. Following the definition of Jy in 
Eqs. (43) and A, in Eqs. (40), we see that if J = 0, then 
Eq. (68) leads to K = 0. Therefore we can interpret 
the parameter J in action (3) as the deviation from the 
standard Ising model in the boundary theory. 

Here some remarks are in order. First, the equality 
Eq. (67) is exact according to the AdS/CFT correspon- 
dence. However, the expressions for the grand thermo- 
dynamic potentials in both sides are just some approxi- 
mations. In gravity side, because we used the probe limit 
which neglects the back reaction of polarization field to 
the background geometry and to external field and we 
only computed the grand thermodynamic potential in 
the classical level which neglects the quantum correction. 
In the boundary side, we used mean field expansion and 
only took the tree level of the quantum fluctuation into 
account. So the relationship (68) is only approximately 
valid and our explanation for J is only a qualitative de- 
scription. Second, the Ising-like model (37) can only de- 
scribe the local moment system, which is usually suitable 
for insulator. For metal magnetic materials, we need to 
use different model. So what the meaning of J is in these 
materials needs to be considered in the future. 

With the grand thermodynamic potential in Eq. (66), 
we can obtain the expression of magnetic moment in the 
ferromagnetic phase as, 


2 
Np? = |D 2 _ — ry) ¥?. (69) 
—16JS pa, 


This confirmed the critical behavior obtained in the nu- 
merical calculations and the critical exponent 1/2 is an 
exact result. We can compute all the coefficients appear- 
ing in (69) and compare them with the numerical ones. 
We can also compute the grand thermodynamic potential 
when magnetic field is nonzero and get the magnetic sus- 
ceptibility and hysteresis loop. The details will be given 
in Appendix B. 


V. PROBE LIMIT BY NEGLECTING THE 
BACK REACTION OF ALL MATTER FIELDS 


A. Spontaneous magnetization and susceptibility 


In the probe limit by neglecting the back reaction of 
all matter fields including the Maxwell field, the back- 
ground geometry is just the planar AdS Schwarzschild 
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black brane with, 


a 3rn 


f(r) =1-3, m (70) 
In the background, the equations for Maxwell field read, 
AJ p? 

2 ee ee 
( r4 ) @ 0, 
2 Po P 
" / 2 
= 0. 1 
o"+=4'- (E+ 2) =0 (71) 


We can directly obtain (r) = f° dr(m? — 4Jp?/r*)p 
with (ra) = 0 and put it into the equation of ¢”, then 
the equation for p can be obtained. So in this probe limit, 
we need only solve the three equations of matter fields 
numerically, 


2 32J pp! 
f 
=0. 
pa (; i A2r4 + 16Jp? — 4m?r4 4 


Note that in this case, there does not exist the AdS2 ge- 
ometry, the near horizon geometry of an extremal black 
brane with vanishing temperature. To have the sponta- 
neous magnetization when the temperature is lowered, 
the restrictions (23) need to be reconsidered. To find the 
restriction about the parameters, let us consider Eqs. (72) 
in the high temperature region where p vanishes and the 
solutions for ¢ and p are the same expressions shown in 
Eqs. (16). We can read off the effective mass square of p 
at the horizon as, 


AT? 9p? 

2 m2 2, 2, _ 2 | 

Meg =M +4Jp(rr) = m 4 ee =m 4 MIAT 
(73) 


Because of J < 0, the temperature term contributes a 
negative term into the effective mass square, which is di- 
vergent when T — 0. Thus we see that in the grand 
canonical ensemble, the instability always appears pro- 
vided that the temperature is low enough. As a result, 
in this case, we need not the restrictions in (23). But, 
the parameters have to satisfy the condition \? < 4m?. 
To see this, let us recall the charge density in Eq. (18). 
In this probe limit, we have, 


A? AnT py 
o(T,u) = 21-7 a) (74) 
Then we reach, 
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We see that it is positive only when A? < 4m?. 

Now we consider the spontaneous magnetization in this 
probe limit. First, let us compute the critical tempera- 
ture Te when B = 0. Similar to the case in the first kind 


of probe limit, the polarization field p is a small quantity 
near the critical temperature, we can neglect the nonlin- 
ear terms of p in the equations of ¢ and p. Then we can 
get (r) = u(1 — ra/r) and p(r) = pr;,/r?m? which are 
identical with expressions in Eqs. (16). The equation of 
p is just Eq. (21). At the horizon, the initial conditions 
are, 


When we perform the numerical computation, we can 
first fix the horizon radius r, = 1 so the temperature is 
also fixed. By adjusting the chemical potential, we shoot 
the boundary condition p, = 0. In the end, we can 
use the scaling transformations (24) and (26) to trans- 
form our results into the case in grand canonical ensem- 
ble where the chemical potential is fixed. As a typical 
example, we also choose parameters as m? = —J = 1/8 
and A = 1/2. The critical temperature is T,./p ~ 1.7871. 
Similarly, we can plot the relationship between p} and 
shooting parameter u, in order to examine whether p 
gets spontaneous condensation when T < T,. We find 
that the solution of source free always appears, which 
results in the spontaneous magnetization of the system, 
and breaks the time reversal symmetry in low tempera- 
tures. 


When the temperature is lower than the critical one 
T, we have to solve Eq. (72) to get the solution of the 
order parameter p, and then compute the value of mag- 
netic moment N, which is also defined by Eq. (14) with 
a = 0. In the top panel in Fig. 5, we plot the value 
of magnetic moment N as a function of temperature. 
We see that when the temperature is lower than Te, the 
non-trivial solution p Æ 0 and spontaneous magnetic mo- 
ment appears. It corresponds a time reversal symmetry 
breaking spontaneously. In addition, let us stress here 
that if the boundary spatial dimension is three, the spa- 
tial rotational symmetry is also broken spontaneously, 
since a nonvanishing magnetic moment chooses a direc- 
tion as a special. The numerical results show that this 
phase transition is a second order one with the behavior 
N x \/1—T/T, near the critical temperature. The re- 
sult is still consistent with the mean field theory descrip- 
tion of the paramagnetism/ferromagnetism phase transi- 
tion. 


Next we calculate the static magnetic susceptibility in 
this probe limit, defined by Eq. (33). Based on the previ- 
ous analysis, the magnetic susceptibility is still obtained 
by solving Eq. (34), and the only difference is the form 
of f(r). Thus we can also set the magnetic field B = 1 
and get m = 5.7(T/Te — 1), which satisfies the Curie- 

xy fi 
Weiss law of ferromagnetism in the region of T > T. Its 
inverse is shown in the bottom panel of Fig. 5. 


12 


03 04 05 06 07 08 09 10 


© Numerical results|] 
— Fitting curve 


1.005 1.01 1.015 1.02 1.025 1.03 1.035 
T/T 


£ 


FIG. 5. Top panel: The magnetic moment N as a function of 
temperature. Bottom panel: The behavior of the inverse sus- 
ceptibility density in the paramagnetic phase near the critical 
temperature. Here, we choose parameters as m? = —J = 1/8. 
The temperature T./: ~ 1.7871. 


B. DC conductivity in the ferromagnetic phase 


The electric transport is also an important property in 
the materials involving spontaneous magnetization. Now 
let us study how the DC conductivity is influenced by 
spontaneous magnetization in this model. In order to 
simplify our computation in technology, we will work in 
the probe limit by neglecting back reactions of all the 
matter fields. This limit can give out the main features 
near the critical temperature. However, in the case of 
near zero temperature, we have to consider the model 
with full back reaction. We will consider this in the fu- 
ture. 

To compute the conductivity, we have to consider some 
perturbations for gauge field with harmonically time 
varying electric field. Due to the planar symmetry at 
the boundary, the conductivity is isotropic. Thus for 
simplicity, we just compute the conductivity along the 
x-direction. According to the dictionary of AdS/CFT, 
we consider the perturbation 6A, = ea,(r)e~***. In the 
probe limit, this perturbation will also lead to the per- 
turbations of polarization field in the first order of e. As 
a result, we have to consider the perturbations for all the 
components of gauge field and polarization field. How- 
ever, if we only care the conductivity in the low frequency 


limit, i.e., T >> w — 0, the problem can be simplified. In 
the low frequency limit, we only need turn on the three 
perturbations, 


OAy = Eaz (je ™, 
Me2=Canie™, (77) 
My = eCty (re, 


and corresponding equations for the three perturbations 
in the low frequency limit read 


mCry AS ppCrx 


Cty — ap Zz + Ow) =0, (78a) 
a, AI ppCry | 
Cra — 3 — Fm? tO(w)=0, (78b) 
2 , 2 , azu? — 
[rf (an — A Cra/4)] + +O(w)=0, (78c) 


r2 f 
with p and p determined by Eqs. (72). Here O(w) is the 
terms with order of w which can be neglected when w > 
0. In general, the term w?/r? f(r) can not be neglected 
since f(r) is zero at the horizon, which leads to that the 
limit of w — 0 is ambiguous. However, at the horizon, if 
we impose the ingoing conditions for Cre, Cty and az, 


Ca= Cc? + Cr =h) +, 
Cra = a [c® a CY) (r = Th) apos iF (79) 
agr =e [a 4 a (r — ra) te] 


with rą = fdr/(r?f), we find the system has a well- 
defined limit when w > 0 if T # 0. At the AdS bound- 
ary with the source free condition, we have the following 
asymptotic solutions, 


Coy = Cry ga E ee, 
80 
On = — ft as- (80) 


2 z tita Ae = Ag+ + — 
rem r 


Then the gauge/gravity duality implies that electric cur- 
rent (Jz) = a,— and the DC conductivity is given by, 


Geo 


lim - . 
w0 War 


c= (81) 

As a holographic application of the membrane 
paradigm of black holes, we can direct obtain the DC 
conductivity from Eqs. (78) using the method proposed 
by Iqbal and Liu in [34]. In fact, the transport coefficients 
in the dual field theory can be obtained from the horizon 
geometry of the dual gravity in the low frequency limit. 
Applying this into U(1) gauge field, this conclusion im- 
plies that the DC conductivity is given by the coefficient 
of the gauge field kinetic term evaluated at the horizon. 
To see this, we assume that T > 0 and w — 0, then we 
can neglect all the terms of w in Eqs. (78). We first note 
that, 


(Jy) A? (Se 
Jim, rf (r)(a%, = Cre /4) = ( v er | 


= —(1— d7/4m?)(J,). 
(82) 
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Eq. (78c) shows that this quantity is conserved along the 
direction r. So at the horizon, using Eqs. (78b) and (78c), 
we have, 


= (1 ~~ A? /4m?) (Je) = lim r’f (dl, a A?’ Cre /4)| 


x AJ ppC;, 
_ 2 _ ; y 
=s |(1 Aa) e+ i. 
(83) 


m? f 
Combining Eqs. (78a) and (78b) and considering the fact 
that Ci, is regular at the horizon, we have, 


16J? p?p? 4J 
fet ER [Ow a A 
MAT m 


T=rp,? 


at r rj. Thus we have from Eqs. (84) and (83) that 
— (1-7/4?) (Je) 
2m? (85) 
m4 + 16J?p?p? /r4) 


+ 
TF 


= lim r’ fal, i rT 


Now let us take the ingoing condition for a, at the hori- 
zon, which tells us that, 


d 


r?’ fal, = -— ar =—iwa,, atr—rf, (86) 
dry 
finally we get, 
iwar(ra) Am? 
Jn) = s (87 
We) = oe, | at arp EP 


Here po and po are the initial values of p(r) and p(r) at 
the horizon, which can be computed from Eqs. (72). In 
the low frequency limit, Eqs. (78) imply that the electric 
field is constant, i.e., lim;=,, @2(r) = dz+. So we have 
seen that we can obtain the DC conductivity near the 
horizon as, 


7 1 Mm? 
~ 1—d2/4m? A(m4 + 16.J?p2 62 /r*) 


o (88) 
In the top panel of Fig. 6, we plot the numerical result on 
the DC resistivity 1/o as a function of temperature. We 
see that the resistivity shows a metallic behavior when 
the temperature is below the Curie temperature. 

With the appearance of ferromagnetism, DC resistiv- 
ity decreases when the sample gets cooling, which shows 
in many interesting phenomena in condensed matter 
physics, especially in a class of manganese oxides which 
are widespread because of the discovery of colossal mag- 
netoresistance (CMR) and receive a lot of interest both 
in theory and experiment [36, 37]. Note that this effect 
has a completely different physical origin from the “gi- 
ant” magnetoresistance observed in layered and clustered 
compounds. In recent twenty years, CMR is among the 
main topics of study within the area of strongly corre- 
lated electron systems and its popularity is reaching the 
level comparable to that of the high-temperature super- 
conducting cuprates. In the bottom panel of Fig. 6, we 
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FIG. 6. Top panel: DC resistivity vs temperature in our 
model. Here we choose parameters as m? = —J = 1/8 and 
A = 1/2. The critical temperature Te/u ~ 1.7871. Bottom 
panel: Temperature dependence of resistivity for various sin- 
gle crystals of Lay_,SrzMnO3. Arrows indicate the Curie 
temperature. For more details, see Ref. [35]. 


show the experimental data from a typical CMR ma- 
terial La,_,Sr,;MnO3 as an example. We see that our 
model gives a very similar behavior as the latter with 
x > 0.175. Of course, we should mention here that there 
still exist some differences between our model result and 
experimental data on CMR. In general, when T > Te, 
the material shows a semiconductor or insulator behavior 
and the DC resistivity increases with cooling the sample. 
In our model, however, the DC resistivity is a constant 
when T > T.. So this model only gives partial property 
of CMR when T < T,. But this is an exciting and en- 
lightening result, because it implies that this model can 
lead to a possibility to build a holographic CMR model 
and to investigate this typical and important strong cor- 
related electrons system in the AdS/CFT setup. We are 
going to investigate this issue in the future. 
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VI. PHASE TRANSITION WITH BACK 
REACTION 


In previous sections, we have studied the spontaneous 
magnetization in two kinds of probe limit. However, 
probe limit may lead to some information lost. For ex- 
ample, in some holographic superconductor models, it 
will lead to the appearance of the first or zeroth order 
phase transition when the strength of back reaction gets 
beyond some values (see Refs. [26, 38, 39], for example). 
In addition, when the temperature is low enough, the 
probe limit may lose its validness. To have a complete 
phase diagram for the holographic model, we need to go 
beyond the probe approximation and to include the back 
reaction. 


A. On-shell free energy 


The model admits various solutions, in order to de- 
termine which phase is thermodynamically favored, we 
should calculate the free energy of the system for both 
normal phase and condensed phase. In gauge/gravity du- 
ality the grand potential Q of the boundary thermal state 
is identified with temperature T times the on-shell bulk 
action with Euclidean signature. Since we are consider- 
ing a stationary problem, the Euclidean action is related 
to the Minkowski one by a minus sign as 


— Sp = J tavar + + Lm), (89) 


where Lm = ~F" Fy + à? Lo and g is the determinant 
of the metric. We first show that, when evaluated on a 
solution, this action reduces to a simple surface term at 
the AdS boundary. From the symmetries of the solution 
(5) and (6), the yy component of the stress energy tensor 
only has a contribution from the terms proportional to 
the metric. Thus, the gravitational filed equations imply 
that 


6 
The Euclidean action is then 


— SE = f arr fer? — A? fpp' e?) 
V2 


B(4B — e 


2 , 


(91) 


r 


The surface term on the horizon vanishes since f(r) = 0. 
Thus we can get the Euclidean action that contains the 
surface term at Too as, 


— Sg 
Və 


= 2r? fer/? = X fop erra 


y2 a/2 
por A pye l 


r2 


As the first item of Eq. (92) diverges when r —> oo and 
must be regulated. This counter terms we need to reg- 
ulate the action are the standard ones (see for exam- 
ple [25]): 


1 
Sua. = gop | Vg 2K +4/L) rars (93) 


where goo is the induced metric on the boundary r = ræ 
and K = g$ V „n, is the trace of the extrinsic curva- 
ture (n” is the outward pointing unit normal vector to 
the boundary). The summation Sensbell = Sp + Set. 
is now finite in the limit r — oo. The regularized ac- 
tion becomes, after considering the asymptotical forms 
in (13), 
Q= T Son-shell /V2 


Euclidean 


®  B(4B — d?p)e%/? 
=- f Ma Ae 4 2fo. (94) 


B. Tensor hairy solutions and phase transition 


We are interesting in the black brane solutions with 
nontrivial space-space component p of the tensor field. 
For this, we have to adopt the numerical method to find 
such solutions. Without loss of generality, the location 
of rz, can be fixed to be unity in our numerical calcula- 
tion. We are then left with two independent parameters 
{p(rn), p(rn)}. By choosing p(r;,) as the shooting param- 
eter to match the source free condition at r > ov, i.e., 
p+ = 0, we finally have a one-parameter family of solu- 
tions labeled by the value of p(r;,). Other coefficients can 
be expressed in terms of those parameters. After solving 
the set of equations, we can calculate the spontaneous 
magnetization N and free energy density. 

With a fixed m?, we scan a wide range of J and À 
within the limitation Eq.(23) in 3-dimension plane in or- 
der to trace out the evolution of critical temperature Te 
versus these parameters. Fig. 7 plots the critical tem- 
perature Te as a function of J and A in the cases with 
m = 1/4, 1/8. 1/16 and 1/160, respectively. We can see 
from Fig. 8 that the critical temperature is weakly de- 
pendent of the parameter A. Moreover, the smaller the 
value of m?, the larger the phase transition temperature 
T, when the same value of J and A are considered. For 
each value of À, the analytical black brane solution (16) 
always exists. However, for sufficiently low temperature, 
we always find additional solutions with non-vanishing 
p which are thermodynamically favored. That is to say, 
for each value of A we take, there is a phase transition 
occurring at a certain temperature Te where the black 
brane developing a new “tensor hair” with nontrival p 
becomes thermodynamically favored. In the dual field 
theory side, it means that magnetic moment acquires 
a vacuum expectation value breaking the time reversal 
symmetry spontaneously. 

The top panel of Fig. 9 presents the condensate as a 
function of temperature for A = 1/8,1/2 and 7/10, from 
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FIG. 7. The critical temperature Te with respect to the model 
parameters J and à. Here (a), (b), (c) and (d) correspond to 
m? = 1/4,1/8,1/16, and 1/160, respectively. In (a), some 
part of the surface is cut because that J and A should satisfy 
the relation (23). 


1.77% 


FIG. 8. The critical temperature T. as a function of » in the 
case with model parameter m? = —J = 1/8. 


which one can see that N rises continuously from zero 
at Te. For small À, the curve is similar to the case with 
the probe limit (compare the case of A = 1/8 in Fig. 9 
and the plot in the bottom panel of Fig. 2), We can see 
from the figure that when A gets larger, the condensate 
increases. However, near the critical temperature, the 
square root behavior still holds as 


N x yT -T.. (95) 


The free energy difference of the condensed phase and 
the normal phase is expressed by AQ, which is plotted 
in the bottom panel of Fig. 9. It is obvious that below 
the critical temperature T., the state with non-vanishing 
magnetic “tensor hair” is indeed thermodynamically fa- 
vored over the normal phase because AQ always be less 
than zero. Moreover, our numerical calculation indicates 
that the order of the phase transition is only second or- 
der, no matter how we increase the strength of the back 
reaction in parameters allowed region, i.e., A?/4m? < 1. 
The first order or even the zero order phase transition 
does not appear in this model. 
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FIG. 9. The magnetic moment N (top panel) and the grand 
potential deference AQ (bottom panel) as a function of tem- 
perature with different A. Here m? = —J = 1/8 for both 
figures. 


VII. SUMMARY AND DISCUSSION 


In this paper we have presented a holographic model to 
realize the paramagnetism/ferromagnetism phase transi- 
tion in AdS black brane background by introducing a 
massive 2-form field in the bulk. This 2-form field cou- 
ples to the background Maxwell field strength and carries 
self interaction. 

The model admits a new analytical black brane so- 
lution with a non-trivial time-space component of the 
tensor field. The properties of the black brane solution 
depend on the value of A?/4m?. When \?/4m? < 1, 
this new solution is very similar to the planar AdS RN 
black hole. But if \?/4m? > 1, the black brane solution 
is chemically unstable in grand canonical ensemble un- 
less the temperature and chemical potential satisfy some 
additional conditions. In that case, there is not corre- 
sponding extreme black hole, i.e., the horizon radius is 
zero when temperature goes to zero. A very special case 
is that when \?/4m? = 1, the space-time geometry is 
just the planar AdS Schwarzschild geometry but both 
the U(1) field and 2-form field nonvanish. 

For our goal in this paper we focused on the case 
with \?/4m? < 1 so that an asymptotic AdS geometry 
near the horizon emerges when the temperature tends to 
zero. By this emergent AdS2 geometry, we obtained the 
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condition in which the spontaneously symmetry break- 
ing can happen. If the parameters satisfy the condition, 
time reversal symmetry will be broken spontaneously and 
the paramagnetism/ferromagnetism phase transition can 
happen when the temperature is lower than a critical 
value. 

In order to understand the properties of this holo- 
graphic ferromagnetic phase transition, we investigated 
the paramagnetism/ferromagnetism phase transition in 
two kinds of probe limit and in the case with full back 
reaction, respectively. 

In the case of the first kind of probe limit where the 
model parameters A — 0. This probe limit neglects the 
back reactions of the 2-form field to the Maxwell field and 
background geometry and makes it simple to study the 
behavior of spontaneous magnetization. In this probe 
limit, we computed the critical exponents by both nu- 
merical and analytical approaches, which are agreement 
with the mean field results. In addition, we obtained a 
Ginzburg-Landau-like free energy near the critical tem- 
perature for the holographic model and argued that it is 
agreement with results from the Ising-like Hamiltonian 
in film system of the condensed matter physics and re- 
lated the model parameter J to the parameters in the 
Ising-like model. 

In the second kind of probe limit, we neglected all back 
reaction of matter fields to the background geometry but 
considered the interaction between the tensor field and 
the Maxwell field. In this probe limit, we are able to 
study the influence of spontaneous magnetization on the 
electric transport properties in a relatively simple way, 
where the background geometry is fixed. We found that 
the critical exponents are the same as ones in the first 
kind of probe limit. We also computed the DC resistivity 
in this probe limit. It was found that the DC resistivity 
is suppressed by spontaneous magnetization and shows 
a metallic behavior. This is very similar to the strong 
correlated phenomenon named CMR effect found in the 
some manganites. 

Next we considered the case with full back reaction 
and solved the full equations of motion numerically. It 
was found that the free energy difference between the 
condense phase and normal phase, AQ, is zero at criti- 
cal temperature and always negative when T < T,. The 
phase transition is always a second order one as one in- 
creases the strength of the back reaction. 

Main calculations are made in 3+ 1 dimensions in this 
paper, but it can be easily extended into the higher di- 
mension case. In the latter case, the space rotation sym- 
metry will be broken spontaneously when the phase tran- 
sition happens. In addition, we can also generalize this 
model to the case with the Lifshitz symmetry in the bulk 
and study the influences of the Lifshitz dynamical expo- 
nent z on the condensate both in the probe limit and in 
the case with the back reaction. In addition, in all the 
calculations in this paper, we assumed that the solution 
in the bulk or the phase at the boundary is homogeneous, 
which is a strong assumption. In fact, inhomogeneous 


solution may exist even in a model whose Lagrangian 
has translation symmetry and the boundary conditions 
are homogeneous (see Ref. [17], for example). In many 
materials, the inhomogeneous phase can appear sponta- 
neously in a chemical homogeneous materials. Therefore 
it is of great interest to study whether the inhomogene- 
ity could appear spontaneously in this model. Although 
we focused on the ferromagnetic phase transition only in 
this paper, which has been understood well in condensed 
matter physics, it offers a framework rather than only a 
specific model, which can be regarded as basic starting 
point to understand more complicated phenomenon in- 
volving spontaneous magnetization. As a result, there 
are various prospects to study and we expect more excit- 
ing results could be reported in the future. 
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Appendix A: Ginzburg-Landau theory for Ising 
universality class 


In this appendix, we will compute the grand thermo- 
dynamic potential of the Ising-like Hamiltonian (38). We 
can rewrite it into the form of \¢* theory, 


1 s34 
H= pes [5025 + ms") + an 


7i (A1) 


Turn into the Lagrangian form” and go to the case with 
Euclidian signature, we have the action, 


1 fF 
Se=3 f ae J ax{(0,3)2 + (Pa? m 
0 


(A2) 
with 8 = 1/T, the inverse temperature of the system. We 
see that the mass dimension of coupling constant [A,] = 
1, so effective action (A2) is renormalizable. In order 
to compute the quantum effective potential, we take the 
mean field expansion method. We split as usual the field 
into classical background part and quantum fluctuation 
part, 


(A3) 
The quadratic part of 7 in the action is then controlled 


by a kinetic operator of the form, 


SE 
525 [Sea] = a 


S= Sa +n. 


V? + M? (Sa), (A4) 


T Here we assume the relativistic depression relation. 


17 


where 
z2 
As S4 


M? (Sa) = —m + 2 


(A5) 
In addition, the Euclidian action contains the linear 
term of 7, but this term contributes nothing and can 
be dropped. And, if we neglect the cubic and quartic 
order terms of 7, we can obtain the result in mean field 
approximation by only taking the tree level of 7 into ac- 
count. This is a good approximation if A, and AsScı are 
small®. Then it is easy to find the grand thermodynamic 
potential density in the mean field approximation as, 


i aa lioa dk w i Bw 
Q= —5Ms Se + gaia +J L23 tT ln(1 — e7®®)]. 
(A6) 
Here w = yk? + M? (Sa). The first term in the integral 


of (A6) contributes the zero-point energy and can be ne- 
glected. Then the second term in the integral of (A6) 
then can be written as, 


Te f= J22 a? 
=f dg ain(1 — e7 V” t*0)], (A7) 
0 


where 22 = B? M? (Sa). The result of (A7) does not ad- 
mit a compact expression. We can get a series expression 
by high temperature expansion, where x9 is treated as a 
small quantity, 


1 dx xIn(1 — e~V +0) = 5 Cnr”, (A8) 
0 n=0 


where 
co = —C(3), & = yE/2 +. 
Here yg = 0.5772--- 


(A9) 


is the Euler constant and ¢(3) = 


1.202- - - is the value of Riemann-¢ function. In this limit 
we have, 
T yem? 1, As yET As 34 
Og 3 T2 i s 2\>2 sl 
EOT + TEMS] + (E — 282, + SH, 
(A10) 


from which we can obtain the expression in (39). 


Appendix B: Semi-analytic calculations near the 
critical temperature 


1. Spontaneous magnetization 


In this appendix, we will compute the values of N1, ao 
and a, appearing in section IV. With these we can get 


8 In the vicinity of critical temperature, \53,; is always small. So 
in this region, the cubic term can always be neglected. However, 
As in general may not be a small quantity. In this case, loop- 
corrections must be considered. 
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FIG. 10. The value of A1 with respective to temperature in the 
case of m? = —J = 1/8 and w = 24 (left) and w = z? (right). 


the coefficients in (69) and compare with the numerical 
results in the previous sections. 

Let us first compute N; and a,. For this we have to 
first find the eigenfunction p1, which is the solution of, 


d 


dpn| _ 
SoZ] +a@m=0 æ 
at T = T, with the conditions, 
pı(1)=1, pip =0. (B2) 


For convenience, here we do not assume that {pn} form 
an unit base. Thus we have, 


1 . 1 i 8 4 
Ni = — d =z d B3 
1 Ci f piaz, ay ae Z paz, ( ) 


here C4 is the normalization coefficient and 


1 
C? = (p1, 1) = wpidz. (B4) 
0 


In order to compute ao, we need to solve equation (B1) 
in the limit T — T; . To clarify that the results are 
independent of the specific form of weight function, we 
choose k = 3,4 as two examples. Then we fit the relation 
At = ao(T/T. — 1) to find ap. Figure 10 shows that 
A; and T/T. — 1 indeed satisfies a linear relation very 
well. Numerical results show that Ni œ 2.0412, a1 ~ 
0.7276,a9 ~ 1.1309 for k = 4 and Nı œ 1.8370, a1 % 
0.4772,a9 ~ 0.9155 for k = 3. We have, 


? 


N?ao 


N?/ u = —= 
—16Jsa u2 


(1 - T/T.) ~ a(1 — T/T). (B5) 


with ag ~ 4.966 for k = 4 and ag ~ 4.964 for k = 3. We 
see that different weight functions give different values for 
Nı, a, and ag, but the same value of magnetic moment 
N(up to a numerical error). 

The value of ap can also be obtained directly by solving 
ODE (52). In the region near the critical temperature, we 
assum A, = aọo(T/Te— 1). Note that all quantities in (52) 
are the functions of temperature, so take derivative with 
respect to T and evaluate at T = T., we get, 


dP adp _ ao 
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Here p is the eigenfunction of (B1). Now treat pr = apy 
as an unknown function to be solved, then the task to 
find a) becomes to solve a non-homogenous eigenvalue 
problem, 


ao dP 


Pos |= B 
Pr=| p gp | ht (B7) 


At the AdS boundary, pr has the same asymptotic be- 
havior as (53), thus we can impose the boundary condi- 
tions as 


lpr(1)| < 00, pr+ =0. (B8) 


We find that pr € D. We then use the basis {pn} to 
expand pr, i.e., 
i dn 
m= fg e (B9) 
n=1 ~” 


Here Cn are the modules of pn. Using the fact A, = 0 at 
T = T, and 


(CI py, Por) = X dn(CT p1, Cn Pon) 
Ti (B10) 
= X dnànôın = dià = 0. 


n=1 


we have, 
( ag dP ) agC? i dP d 
=- — = — wpı —pıdz 
P1, T. dT PA T, i Pl dr P1 
apC? 1 dL 
— Pl pidz = 0. 
Te o dT 
(B11) 
Furthermore we get, 
Te f, d 
=~ dzp,— p1. B12 
ao C? I ZP1 ar’! ( ) 


It is very useful to find its equivalent form in the case 
with fixed r, = 1, since it is convenient when we perform 
numerical computation. If we fix r, = 1, the shooting 
parameter is chemical potential u. The relation between 
temperature in grand canonical ensemble and chemical 
potential is given by, 
2 
pat (B13) 
Ary 


Thus the expression (B12) can be rewritten as 
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Te dp? f' dL 
ao = H J dzpı duz”! 
0 H 


uC? dT ras 
8usnT. ft dL 
"307 J, OP gua? 
1 P= [Me 


Here pe is the critical chemical potential when we fix 


a ee Saat . 
rr, = l and T, = rae is the critical temperature in 


grand canonical ensemble. Combining (B3) and (B14), 
we have, 


1 1 L 
Tiel elf pidz)? i dzpı TPL, 


N? / 2y = = 
/ 24 Jp i. dzz4 pt 


1—T/T.). 

(B15) 
We can see that it is determined by the equation (B1), 
but independent of the weight function! The expres- 
sion (B14) depends on the weight function, because it 
depends on C4. In order to check the formula (B14), let 
us compute the values of ap in the cases of k = 4 and 
k = 3. The results gives 1.1286 and 0.9140, respectively. 
We see that, up to numerical errors, they are the same as 
what we have obtained by fitting the curve in figure 10. 
In addition, the expression (B15) gives ag = 4.9560, 
which is very close to what we obtained in the numerical 
calculation. 


2. Susceptibility and hysteresis loop 


When B Æ 0, the susceptibility for T > Te is defined 
as, 


ON 
= li lg Bl 
X= Bo (S5). oo 
In the case with T > Tẹ and B — 0, we can neglect 


the non-linear term, i.e., setting J ş = 0. The solution of 
equation (43) can be expressed as, 


~ B 
p= > Sa m2” 


n=1 


(B17) 


Taking into account the equation (43) with J, $ = 0, we 
have 


a >. 4BJu?z* 
O=Lp+ B=) alp E 
~ a 4BJp?2^ aes 
— zZ 
— Ss" aC; Th wpr = 2 G 
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Multiplying a factor pn/Cn and integrating the above 
equation from 0 to 1, we can obtain 
1 4Ju?z* 


Boyn 
An? 


Thus we can get the magnetic moment density as 


= with Yn = dz. (B19) 


B ty Ng 
N/M = B 5 B2 
[A 2m? = Wn? eed) 
and the magnetic susceptibility 
1 L mN, 
2 niin 
- : B21 
X/A 2m2 = 2Xn ( ) 
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When T > T+, we have Ay = ao(T/T.—1) > 0*. Thus x 
is dominated by the first term in the summation of (B21) 
and its inverse can be expressed as 


2ao 


1), as T =T. (B22 
P ) (B22) 


APX he = 


(T/Te 


In the case of m? = —J = 1/8, we have A?y~1/pe ~ 
4.0520(T/Te — 1), which is very close to our numerical 
result \?x~1/pe ~ 4.0499(T/T. — 1) given in the numer- 
ical calculation. 

Next let us move to the case with B Æ 0. In this case 
from (43) we have, 


1 
f Pn(wPp + B — J;p?z*)dz = 0. (B23) 
0 

Consider (60), we can rewrite it to 

— q2),_ 7 
i pn(wPp+ Bil — me —Jyp°z®)dz =0. (B24) 

0 
Using the expansion expression (60), we have, 
1 
CnC? An + Bryn — | PnJ 5 p2zedz =0, n=1,2,---. 
0 

(B25) 


For convenience, we assume that {pn} is an unit base, 
i.e., C,=1. The equations (B25) are equivalent to (43) if 
we take all the terms in (60) into account. In the case of 
T — T; , assume that the first term in (60) dominates 


only, i.e., |c1| > cn for n > 2 in (61), we have 


B 
N/ = zya ~ &N1/2. 


Taking n = 1 in (B25), we have, 


(B26) 


1 
ah+Bu-eJ; f przdz = oA,+ By —4c} Jfar =0. 


$ (B27) 
For a given temperature T —> Te, we can combine (B26) 
and (B27) to obtain a relation between external magnetic 
field B and magnetic moment density N. Figure 11 shows 
the results with T = 1.057.,T = 0.97, and T = Te, 
respectively, in the case of m? = —J = 1/8. We see 
that it is very similar to what we have obtained in the 
previous work [23]. 
Finally, let us notice that in GL theory, the equation 
for magnetic moment density is, 


A(T —T.)N + A&N’ — B=0 (B28) 


with two positive coefficients A, and Aj. However, it 
is easy to see that the equation for N in our model is 
different from the usual form (B28) from the GL theory, 
which can be obtained by combining (B26) and (B27) 
to eliminate cı. Namely, alhought our model gives the 
similar results near the critical temperature in GL the- 
ory, there exist some differences between the holographic 
model and the GL theory even in the region near the 
critical temperature. 


FIG. 11. The relation between magnetic moment density N 
and external magnetic field B in the cases of T = 1.05T.,T = 
0.9T. and T = Te, respectively. 
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